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Abstract 



Since b quark mass is not asymptotically large, chirally enhanced correc- 
tions which arise from twist-3 wave functions may be important in B decays. 
We thus evaluate the hadronic matrix elements with the emitted meson de- 
scribed by leading twist and twist-3 distribution amplitudes After 
summing over the four "vertex correction" diagrams, we obtain the results 
with infrared finiteness which shows that chirally enhanced corrections arise 
from ^p{x) can be consistently included in QCD factorization. We also briefly 
discuss the contributions from "hard spectator" diagrams. 
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It is well known that two-body, non-leptonic charmless B decays are crucial for extracting 
CKM matrix elements. However, due to our ignorance on how to calculate hadronic matrix 
elements, we conventionally resort to Bjorken's color transparency argument which lead 
to " naive factorization assumption" , 

{MiM2\Q\B) = (M2|Ji|0)(Mi|J2|fi), (1) 

This assumption makes the hadronic matrix elements scale- independent. Noting that Wilson 
Coefficients are scheme- and scale-dependent, the theoretical calculations on the branching 
ratios would then be scheme- and scale-dependent which is unacceptable. To save factor- 
ization hypothesis, a phenomenological parameter Neff is introduced Q which is commonly 
called generalized factorization. However this approach is not satisfactory yet because in 
principle N^fj is process dependent. 

Recently, Beneke, Buchalla, Neubert and Sachrajda proposed a promising method: 
in the heavy quark limit, they show that the emitted light meson can be described with lead- 
ing twist-2 distribution amplitude, the infrared divergences of the hard-scattering amplitudes 
are canceled after summing over the four "vertex correction" diagrams (Fig.(a)-(d)), which 
is a one-loop demonstration of Bjorken's color transparency argument [^. In the heavy 
quark limit, they show that the hadronic matrix elements can be expressed as 

(M1M2IQI5) = (M2I Ji|0)(Mi| J2I5) ■ [1 + Sr„< + O(AqcdM)]. (2) 

If power corrections in l/rrif, can be safely neglected, then everything is perfect. At the zero 
order of as, it would come back to "naive factorization", and at the higher order of a^, 
the corrections can be systematically calculated in Perturbative QCD, which means that 
the decay amplitudes of B meson can be computed from first principles, and the necessary 
input are heavy-to-light form factors and light-cone distribution amplitudes. But in the real 
world, bottom quark mass is not asymptotically large (but 4.8 GeV), and numerically power 
suppression may fail in some cases. An obvious and possibly the most important case is 
chirally enhanced power corrections. As pointed out in ref 0, numerically the enhanced 
factor Tv = — — T — 1.18 which makes the power suppression completely fail. This 
parameter is multiplied by and ag, where is very important numerically in penguin- 
dominated B decays. So an evaluation of the hadronic matrix elements including chirally 
enhanced corrections may be phenomeno logically or numerically important. In this letter, 
we shall examine this problem in some detail. 

Chirally enhanced corrections arise from twist-3 light-cone distribution amplitudes, gen- 
erally called $p(x) and ^a{x). For light pseudoscalar mesons, they are defined as P] 

{P{p'My)^l5q{xm = du e'^-^P'-y+^^'-^^Mu), (3) 

Jo 

{Pip'My)cr,.l,qixm = ^U^^,{p'^z^ - p'^z,) du S^v'-y+^v'-^)<tM^ (4) 

where /ip = ^ /m^ ' z = y — x. We notice that in Ref. color transparency is demonstrated 
in one- loop level in the heavy quark limit. If we want to include chirally enhanced corrections 
consistently, we should describe the emitted light meson with leading twist-2 and twist-3 
distribution amplitudes, which means that we should show the infrared finiteness using 
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twist-3 distribution amplitudes after summing over the "vertex correction" diagrams. In 
this paper, we shall restrict ourselves to ^p{x) while postpone the discussion of ^a{x) to ref 
1^ because of the complicated derivation in proving the infrared finiteness of the "vertex 
correction" diagrams using 

We notice that in Ref [§, the authors have used twist-3 distribution amplitude ^p{x) to 
calculate the strong penguin corrections (Fig.(e)-(f)). The difference of our work from ref 
1^ is that we calculate "vertex correction" diagrams and show the infrared finiteness of ag 
and ag at the order of ag- 
in the following, we take leading twist and twist-3 wave functions ^p{x) to describe the 
emitted light meson, and we will show the infrared finiteness under this approach. 
The |A5| = 1 effective Hamiltonian is given by |^ 



Gp 
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^ V, (criMiii) + CMQlili) + E CMQm\ - Vt{Cr,Qj, + CsgQsg)] + h.c, 

lg=u,c \ fc=3 / 

(5) 



where Vq = V^6V^'^(for b ^ d transition) or Vq = V^^Vg* (for b ^ s transition) and Cj(/i) are 
Wilson coefficients which have been evaluated to next-to-leading order approximation. The 
four-quark operators Qi are 

Ql = {Uaba)v-A{.qpUp)v-A Ql = {Caba)v-A{qpCi3)v-A 

Q2 = {UaUa)v-A{qpbf))v-A Q\ = (c^Ca) v_A(g/3&/3) V-A 

Q3 = {qaba)v-A T.{qpqp)v~A Qi = Y.{qpqp)v-A{qJ>a)v-A 

q' q' 

Qb = {qaba)v-AY.{qpqf})v+A Qe = -2J2iq(3qf3)s+piqaba)s-p (6) 

q' ql 

Ql = U^Mv~AEeq'{qpq'p)v+A Qs = -^lleq>{qpqp)s+p{qJ)a)s-P 

Q9 = l{qaba)v-AJ2^q'{qi3qi3)v-A Qio = lJ2^q'{qi3q[^)v-A{qaba)v-A 

q' q' 

and 

<577 = ^mbqa(T'"'{l + 75)baF^u, QsG = ^rrn,q^a>'''tlpbpGl,, {q = d or s). (7) 

The amplitude of the decays of B to two light pseudoscalar mesons in QCD factorization 
can be written as: 

AiB ^ AfiMa) = % E E uX{M,M,\Q,\B)f, (8) 

where Vp is CKM factor and {MiM2\Qi\B)F is the factorized matrix elements. 

We calculate QCD coefficients with the emitted mesons described by light-cone distri- 
bution amplitudes. For instance, let we consider the "vertex correction" diagrams. Twist-3 
distribution amplitude ^p{x) makes no contribution when considering {V — A) ^(V — A) 
and {S + P) ^{S — P) currents because of their Lorentz structure. As to (V — A) 0(V + A) 
current, there is some subtlety in regularizing the infrared divergences. If we use dimension 
regularization, the infrared finiteness would not hold for the case of ^p{x) after summing 
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over those four "vertex correction" diagrams. That is because wave functions are defined in 
4- dimensions, it may be unconsistent to extend its usage directly to D-dimensions. Thus we 
assign a virtual mass to the gluon propagator and regularize the infrared integrals in four 
dimensions. Then the "vertex correction" contributions of {V — A) (^{V + A) current to 
{S + P) 0(5' — P) operator is infrared finite: 

^=ir^ dx^p{x){m\og^ logx+^logx-Li2— +^^2 -3 }, 9 

47rA'Jo XX X X x 

where a; = 1 — x and Li2{x) is dilogarithm function. This means that we can consistently 
include twist-3 distribution amplitude ^p{x) in the framework of QCD factorization. The 
explicit expressions of of for i = 1 to 10 (for symmetric light-cone distribution amplitudes 
of light pseudoscalar mesons) are obtained as: 

2 SO 
HGmAsp) - 3)^1 + {GmM + Gm.(1) - fu, - fj, + y)C3 

-^C4 + (SGmM + Gm,{sc) + G'm,(1))(C4 + Ce) + Gm^sCs^ , (13) 
«^-^^ + |-S^^^-S^{^^^-^^3./(C. + Ce))ln-^ 

H^G'^M + G'm,(«c) + G'^^imC + C) - If, + ^Cscj , (15) 

"« = ^«+7V-4^Ar'^^ 



+(^ - G'„,(Sc) + iG'„,(l))(^ + C5) - |c,,} , (17) 
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iio - <-io + iTf + TZ^^Q^ 



N An N 

+(-Gm.(^c) + %^)(C^5 + ^) - ^C^t.Gm^s} . (19) 

Here A^" = 3 (/ = 5) is the number of color (flavor), Cp — ^2n^ factor of color, 

Sp — nip/ ml for p — u,c and we deflne the symbols in the above expressions as: (most of 
them are as the same as Beneke's except for G'ji^^{s) and G'{s,x)) 

F = _121nA_i8 + /^ +/^ (20) 
nib 

1 1 
/m2 = / dx g{x)^M2{x), Gm2,8 = J dx G8{x)^M2{x), (21) 



1 

Gm2{s)^ I dxG{s,x)^M2{x), (22) 



1 

G'^^{s)= f dxG'{s,x)<l>l,^{x), (23) 



here $M2(a;)($M2(^)) leading twist (twist-3) wave function of the emitted meson M2, and 
the hard-scattering functions are 

1 — 2x 2 
y(a;) = 3--^lnx-3i7r, G^ix) ^ -—, (24) 

1 

G{s,x)^-4:j duu{l-u)ln{s-u{l-u){l-x)-ie), (25) 


G'{s,x) ^^G{s,x). (26) 

As to ///^ which labels the contributions from the hard spectator scattering diagrams 
(Fig.(g)-(h)), We take the wave function of B meson as 75(^5 — Mb)^b{0 flnd that, 
when considering twist-3 ^p{x) distribution contributions, the hard spectator scattering con- 
tributions are proportional to : 

(i) for the case of {V - A) ® {V - A) and {S + P) ® {S - P) currents , 

/ Mi) / / JjjfaW ^ - g - ^)) I (27) 

I i I y i \ ^ x(i-x) J 

(ii) for the case of (V — A) ® (V -\- A) currents. 
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" df / d. ^'^^ / dr .2o^ 

^ J Mb y J x{l - x) ' ^^^^ 

Under the symmetric distributions of final state light mesons, the logarithmically divergent 
integrals are canceled after summing over two hard spectator scattering diagrams. As a 
result there is no hard spectator scattering contributions to and ag, the contribution of 
hard spectator scattering to other af is as the same as ref [p|,|To|: 

,11 _ 47r2 f^jB } } ^^Ahix) } J.. ^^M^iy) 



In ref , the authors have discussed the contributions of asymmetric distributions and find 
that numerically this effect is very small. In our case, when considering the asymmetric 
distributions, there would appear divergent integrals, but in this case the asymmetric dis- 
tribution corrections would also be small if we parametrized the divergent integrals as an 
unknown parameter (as what have done in ref |^) and could be safely neglected. 

We notice that the above approach of evaluating "hard spectator" contribution is naive. 
For instance, the scale of "hard spectator" contribution should be different from the "vertex 
correction" contribution. While it seems reasonable, for the "vertex correction" diagrams, 
to take the scale /i ~ 0{mb) to avoid large logarithm a^log^, a natural choice of the 
scale of "hard spectator" contribution may be around (9(1 GeV) because the average mo- 
mentum square of the exchanged gluon is about 1 GeV^. Another disturbing feature of 
"hard spectator" contribution is that, as have been pointed out in ref when including 
the contribution of there would appear divergent integral Jq dx- even if the symmetric 
distribution amplitude is applied. This divergent integral implies that the dominant contri- 
bution comes from the end-point region, or in another word, it is dominated by soft gluon 
exchange. However the transverse momentum may not be omitted in the end-point region, 
if so, the corresponding divergent integral would then changed to: 

dxd'kT^k^. (30) 

As an illustration, we do not consider the kr dependence of wave functions (though it is 
certainly not a good approximation), then the above integral proportions to: 



/ 



dxdk^ f dxdy ,^ , 

0^ / • (31) 



x^m^ + k^ J X + y 



The above integration is convergence now, furthermore it is not dominated by end-point 
contribution. This illustrates that the treatment of "hard spectator" diagrams may need 
further discussions. 

In summary, we consider some chirally enhanced corrections arise from twist-3 distri- 
bution amplitudes ^p{x). To include chirally enhanced corrections consistently in QCD 
factorization, we describe the emitted light meson with leading twist and twist-3 ^p{x) 
distribution amplitudes and show the infrared finiteness of "vertex correction" diagrams 
(fig.(a)-(d)). We also briefly discuss the disturbing "hard spectator" contributions. 
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TABLES 



QCD 
v^oemcients 


H = 5.0 GeV 




= 2.5 GeV 

ivT n 


T n 




1.043 + 0.012i 
0.034 - 0.075i 


1.019 
0.184 


1.067 + 0.025i 
-0.011 -0.104i 


1.039 
0.104 


03 
a| 
al 
as 


0.007 + 0.002i 
-0.030 - O.OMi 

-0.035 - O.OOGi 
-0.007 - 0.003i 
-0.047 - 0.003i 
-0.049 - 0.006i 


0.003 
-0.029 

-0.029 
-0.005 
-0.044 
-0.044 


0.010 + 0.004Z 
-0.033 - O.OlSi 

-0.040 - O.OOGi 
-0.009 - O.OOGi 
-0.050 - 0.003i 
-0.053 - O.OOGi 


0.005 
-0.041 
-0.041 
-0.010 
-0.052 
-0.052 


07 X 10^ 
r^a^ X 10^ 

r^a^ X 10^ 
ag X 10^ 
afo X 10^ 
afo X 10^ 


ll.l + 2.7i 

49.0 - 3.0i 

47.3 - 4.6i 
-939.8 - 13.3i 

22.4 + 57.8i 

19.1 + 63.0i 


9.1 
45.3 
45.3 
-913.3 
-135.7 
-135.7 


2.8 + 5.5i 
57.0 - 1.2i 

5G.3 - 1.9i 
-9G9.6 - 25.4i 
G8.1 + 88.9i 
66.2 + 91.9i 


3.6 
49.2 
49.2 
-941.4 
-68.0 
-68.0 



TABLE I. The QCD coefficients a^(7r7r) at NLO and LO for the renormahzation scales at 
H = 5 GeV and n = 2.5 GeV 
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FIGURES 




(e) (0 (g) (h) 



FIG. 1. Order of corrections to hard-scattering kernels. The upward quark Hnes represent 
the ejected quark pairs from b quark weak decays. 
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